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PACS 78 . 67 . W j - Optical properties of graphene 

PACS 73 . 20 . Mf - Collective excitations (including excitons, polarons, plasmons and other charge- 
density excitations) 
PACS 71.35. Ji - Excitons in magnetic fields; magnetoexcitons 

Abstract. - We show that graphene in a strong magnetic field with partially filled Landau lev- 
els sustains charged collective excitations — bound states of a neutral magnetoplasmon and free 
particles. In the limit of low density of excess charges, these are bound three-particle complexes. 
Some of these states are optically bright and may be detected in spectroscopy experiments, pro- 
viding a direct probe of electron-electron interactions in graphene. The charged excitations can 
be classified using the geometrical symmetries — non-commutative magnetic translations and 
generalized rotations — in addition to the dynamical SU4 symmetry in graphene. From the SU4 
symmetry point of view, such excitations are analogous to bound states of two quarks and one 
antiquark qqq with four flavors. We establish a flavor optical selection rule to identify the bright 
states for experimental studies. 



> 

rn 

in 

o 
o 



X 



■ Since the isolation of graphene, the unusual behavior of 
its electrons in a magnetic field has been the subject of in- 
tense study [TII2]. Early on, the anomalous integer quan- 
tum Hall effect was observed, with plateaus detected at 
filling factors v = ±2(|n| -I- 1), where n is the Landau level 
(jLL) index |21ll]. Subsequently, in higher fields, plateaus 
were seen at i^ = 0, ±1, ±4, signifying a lift in the fourfold 
(^pin and valley) LL degeneracy, most likely due to the 
many-body effects [5] . Recently, fabrication of high qual- 
ity samples has allowed the fractional quantum Hall effect 
(■FQHE) to be observed at z^ = ^, indicating again the 
presence of a strongly correlated electron state [51[7] . Im- 
portantly, such signatures of correlated many-body physics 
emerge already at low magnetic fields, and in higher fields 
these correlations determine the physics. Here we estab- 
lish the existence of collective charged excitations which 
are essentially bound states of a neutral magnetoplasmon 
binding free particles. In the limit of low density of ex- 
cess carriers such collective excitations can be thought of 
as three-particle interacting states. We will demonstrate 
that in graphene complexes of electrons, akin to qqq quark- 
triplets in high-energy physics, arise naturally. This is be- 
cause the spin and valley pseudospin in graphene combine 



to form four flavors analogous to the four flavors of 1st 
and 2nd generation quarks [HIS]. Analysis of the sym- 
metry content of these states allows the classification of 
optically active and dormant states, highlighting the use- 
fulness of high energy analogues in graphene. 

We study the formation of such excitations in pristine 
graphene in strong magnetic fields with a low density of 
excess electrons in partially filled LL's. More specifically, 
we consider inter-LL excitations which, in the presence of 
excess charges, may bind one additional electron (or hole). 
Such excitations (i) lie above charged topological collective 
excitations, Skyrmions [T0lI12) . and (ii) are directly acces- 
sible via optical and tunneling experiments. In particular, 
we are interested in the behavior of graphene close to the 
charge neutrality point and thus focus on filling factors 
V = ii±e for /I = 0, ±1, ±2 and e <C 1. By way of illustra- 
tion, let us consider the filling factor v = —2 -|- e, so that 
the n — Q LL is almost empty (see Fig. [TJ inset). A neu- 
tral collective excitation, or magnetoplasmon, is formed 
when an electron, excited to an empty state in the n = 
LL, becomes bound to the hole formed in the n = — 1 LL. 
On their own, such neutral excitations form a continuum 
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Fig. 1: (color online) Predicted absorption spectra in pris- 
tine graphene with v = — 2 + e for the bright state at energy 
~0.09£'o. The inset schematically shows the relevant optical 
transitions. Energies are given in units of Eq (see text) with 
zero energy corresponding to the CR energy, hojc- 



of extended states, with the lower continuum edge corre- 
sponding to the cyclotron resonance (CR) mode at energy 

huc — vf \ ^'^^^ [T5] . In addition, they may bind one of 
the few excess electrons in the n = LL, thus forming a 
negatively charged trion (ee/i) complex, X~ . These dis- 
crete bound states are stable when located below the mag- 
netoplasmon continuum. We demonstrate the existence of 
such states, and determine their spectral properties. We 
find that the symmetry of graphene severely restricts opti- 
cal transitions so that most collective states are dark. The 
bright bound states give rise to additional optical absorp- 
tion peaks below the CR as shown in Fig. [T] and whose 
observation is within experimental reach. The binding en- 
ergies of the extra electron in the X~ are ~ O.lE^Oj similar 
to large FQHE gaps expected in pristine graphene P^fTS] 
(^0 - 40meV at B = 10 T). Here Ea = y^^ is the 
characteristic energy of Coulomb interactions in strong 
fields, Eq ^ hujc, Ib = (hc/eB) ' is the magnetic length 
and £ the effective dielectric constant. In this regime, vir- 
tual transitions between LL's are suppressed as powers of 
Eq/HoJc and are neglected, constituting the high magnetic 
field approximation P^[ni [TMT5 ] . 

Due to the particle-hole symmetry in graphene, posi- 
tively charged trions X+, or (ehh) states, exist as X~ 
counterparts. This occurs for the complementary filling 
factors —ly, such as, —v = 2 — e, when the n = LL is 
almost completely filled but contains a few holes. The X+ 
states at — i^ have the same energies and optical strengths 
as their X~ counterparts at v, but are active in the oppo- 
site circular polarization. 

With the underlying dynamical SU4 symmetry of 
graphene, resulting from two possible spin (t, 4,) and two 
valley pseudospin (ft", -(I-) projections, the trion states are 
in fact completely analogous, from the symmetry point of 
view, to bound states of two quarks and one antiquark 



with four flavors, qqq for the X^ and qqq for the X^ 
[5] . The corresponding multiplets are shown in Fig. [5] LL 
splittings in Fig. HJa) are shown for clarity. In particular, 
the state designated in Fig. [T] as X^ realizes the funda- 
mental representation of the SU4 group. Fig. [Ifb). The 
elementary particle physics classification of states and the 
techniques of the SU4 group ^ turn out to be very useful 
in graphene. They allow us to perform a complete group- 
theoretical analysis of the bound states and to determine 
the multiplicities of the emerging groups of excitations. 
Importantly, these enable us to establish a flavor optical 
selection rule, explaining which of the states are dark and 
which arc bright for experimental studies. 

In the following, we focus on the X~ states, which may 
occur for fillings v — fJ-+£, where e ^ 1 and /i is an integer; 
the X~^ case proceeds analogously. The operator creating 
a charged collective excitation is 
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is used for 



all completely filled levels. The composite index N = 
{h^SztCTz} carries the LL index n, the spin Sz — ±^ and 



the pseudospin a^ — ^\ (labeling valleys K and K'), 
while N — {Ni,N2,N3}. The oscillator quantum num- 
ber, m, represents in the symmetric gauge A = ^B x r 
the guiding center for a single particle, yielding the LL 
degeneracy. As for the neutral case [I^, several of the 
-^Ij N states having same k, J^ and equal or close ener- 
gies (but different N's) are strongly mixed by Coulomb 
interactions. As a result, the true excitation is given by a 
linear combination of these (see Fig. [2|). 

The geometrical symmetries of the system allow the clas- 
sification of the composite excitations by two orbital quan- 
tum numbers, half integer Jz and integer fc = 0, 1, . . .. The 
former is due to the symmetry under generalized rotations 
in graphene. It corresponds to the total Jz = ^^jiz, 
where a single-particle operator jiz — hz ~ \'^2 ® tiz in- 
volves the orbital liz and the sublattice isospin Tiz angular 
momentum projections [181I19J . The oscillator quantum 
number, k, is due to magnetic translations whose genera- 
tor K = X^i -^i commutes with the Hamiltonian H . The 
single particle operators are K^ = Tr^ — — f^ x B, where 
TTi is the kinematic momentum operator and e^ = ±e the 
particle's charge. For charged complexes, the components 
do not commute: [KxiKy] — iheB/c for the X~. How- 
ever, using the mutual commutativity of H, Jz, and K^, 
the composite states can be chosen to be eigenstates of 
K2 with eigenvalues (2k + l){h/lBf and fc = 0. 1. . . . pp]. 
Physically, the oscillator number k determines the distance 
at which the guiding center of a charged complex as a 
whole is located relative to the origin. For a uniform sys- 
tem this leads to the Landau degeneracy in k. Therefore, 
charged composites occur in families starting with a seed 
state, which has k = Q and a certain value of J^; the lat- 
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y = -2 + c 





Fig. 2: (color online) (a) Filled/empty circles (•/o) denote electrons/holes; hatched regions, the Dirac sea. States in the boxes 
are resonantly mixed by the Coulomb interactions, (b) Quartet [4] and (c) [20], and (d) [36] SU4 multiplets describe the trion 
X~ states. Shading denotes the different levels of charm. The neutral (eft.) state in (b) Qo = | [dd + uu + ss + ccj is the SU4 
singlet. The electron in the n — 1 LL shown in (a) is resonantly admixed to the X~ from the [4]-multiplet with Jz = |. 



ter cannot be guessed a priori. The offspring k — 1,2,... 
are generated by successive action of the raising ladder 
operator K^ = -JM— 



y I . They have decreasing 



eigenvalues of Jz equal to J^ — k because of the relation 



Jz,K_ 



-K. 



In addition to having the same ener- 
gies as each other, states in a given family also exhibit the 
same optical properties j21j . 



The Hamiltonian is given hy H = Hq + Hi^t with the 
interaction part Hint — H^e + Hhh + H^h and Hq — 
Y.N,m ^NC^N,nCNm-Y._N,m eArd]v„djv„j being the non- 

interacting part. Note that the range over which a given 
N quantum number is summed, depends on whether it 



indexes an electron c'j^.^ or hole djy^ operator. The 
single-particle energy cn = eN + Ese{N) denotes the bare 
LL energy e^v renormalized by e-e exchange self-energy 
corrections £'sE(iV) [T71 UU |22] • The bare LL energies 
e^v = sign(n)?ici;c\/|ri| -I- fkUsSz + fujjyaz may phenomeno- 
logically include the Zeeman splitting fvjjs and a possi- 
ble valley splitting fia;„; these are assumed to be small, 
fiuJs,v ^ Eq. The electron-electron iJee, hole-hole Hhh, 
and electron- hole Hf-h parts of the interaction Hamiltonian 
Hint are found by projecting onto particular LL's. Their 
forms for excitations at 1/ = — 2 + e are (cf. Fig. [H^a)) 



H, 



- V V W^^'"'i^^™^c^ c'^ c c 

2 Z_^ Z_^ '^NimiN2m2^N[m\^N:^'m'^^N2Va2^Nivai 



nil ,"^2 si ,cri 



(2) 
with N, ^ Nl ^ {0,Si,ai} and i = 1,2. The 

graphene Coulomb interaction vertex VV^j^™^' jy^™^ ~ 
^si.s' S(7i. a' Ss2,s' ^<T2.<7' coHserves spin and pseudospin and 
may be expressed in terms of the vertices calculated using 
the standard 2D electron wavefunctions in LL's [IH]. The 
Hamiltonian of the e-h interactions consists of two parts. 



H, 



eh 



m 



HI\ , where 



hl=ee <i:i 



N[m[N^m'-\ ,) , 



7t 



nil ,ni2 si ,(Ti 



N2m2 N[ni{ N^ni'2 N2m2 Ninii 



with Ni ^ N[ ^ {0, si,cri} and N2 
The e-h vertex 

.^N[ni'i N'2m'2 _ y^N'iin'i N2m2 
''NimiN2ni2 ^Ninii N^ni'^ 



(3) 



N^ = {-l,S2,a2}. 



W^r'TZC^ (4) 



Ninii Nl2m'2 

describes the direct e-h Coulomb attraction and the ex- 
change e-ft, repulsion. The term iJ^I^ — SHl\ + H.c. with 



sHii = E E <:; 



nil ,ni2 si ,ai 
m.' .m' «2,CT2 



N2m2 t t jt 

. ^2"!^ '^N[7n['^N^ni'2 N2ni2'^Ninii 



(5) 
and A''i = {l,si,cri}, N[ = {0,si,cri}, N2 = {-1, 32,(72}, 
and N2 — {0, S2,o'2} does not conserve the number of 
particles in each of the LL's involved. It describes the 
resonant conversion between a three-particle trion state 
and a single electron state in the n = 1 LL and con- 
tributes ~ Eq. These indicated interaction terms exhaust 
all Coulomb contributions ^ Eq for the considered charged 
excitations at z^ = —2 -I- e. For other excitations and at 
other filling factors, the effective interaction Hamiltonians 
are constructed analogously. 

The excitations may additionally be characterized by 
their total spin projection Sz and pseudospin projection 
(valley index) E^ [I1I2- However, there exists a larger 
SU4 group, which includes the direct product of the spin 
and pseudospin groups SU2 (^ SU2 C SU4 as a subgroup. 
More precisely, the Hamiltonian of the Coulomb (and in 
fact any long-range) inter-particle interactions is SU4 sym- 
metric up to small symmetry-breaking terms |16j . The 
dynamical SU4 symmetry in graphene is based therefore 
on the equivalence of the two valleys (tl,-!!) and the two 
spin states (t)4-)) including the continuous interchange be- 
tween spins and valleys. The four basis states (flavors) 
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that form a fundamental representation of the SU4 group 
are 4.JJ., li[, tJJ-, and ffT- We shall label them in anal- 
ogy with the four flavors of down, up, strange, and charm 
quarks, by {i-D-, ift, t-D-, ttT} = {d,u,s,c}. The genera- 
tors of the SU4 group can be constructed using the bi- 
linear combinations conserving the number of electrons 
and holes. We have Cij = J^Af^^iCA/j - Z^aA ^Jvj ^A^* ■ 
where the composite index Af = {n,m} and i,j G 
{d,u, s,c}. The generators satisfy the commutation re- 
lations [C'ij, Cki] = SjkCii — SiiCkj- The operators of total 
spin S and pseudospin S can be expressed in terms of 
the Qj. For example, J^^ ^ ^ {Ccc + Cuu - C'ss - Cdd) 
and E+ = E^^ -I- iSj, = Cud + Ccs- The degeneracy of the 
eigenstates follows from the invariance of the Hamiltonian 

H = WHU, where U = exp (i J^ij ©y C'y ) , with 8^ the 

SU4 transformation "angles" [9l[T6]. 

To determine the multiplicities of the 64 possible eeh 
flavor states, we need to decompose the direct product of 
the three SU4 multiplets. The electron states form the fun- 
damental SU4 [4]-multiplet with four states represented 
by the Young diagram Q. The hole (antiparticle) trans- 
forms as a conjugate multiplet [4] and is represented by 
the Young diagram H, which is one box short of the uni- 
tary singlet -. The latter transforms as the SU4 vacuum. 

The decomposition of the direct product of the SU4 eeh 
states gives [4] ® [4] ® [4] = [4] © [4] © [20] © [36]. The 
corresponding Young diagrams are, respectively. 



bound trions X are shown for filling factors v = ±2 -I- e 
in Fig. [21 We truncate the basis in Eq. ([T]) , so that only 



□0n®_ 



(6) 



Graphically, the resulting multiplets with their flavor con- 
tents resolved are shown in Fig. [2rb)-(d). There, we use 
the standard elementary particle physics SU4 representa- 
tion [nilUI- Namely, the three orthogonal axes show: the 
flavor isospin projection T^ — ■^ {Cuu ~ Cdd), the hyper- 
charge Y ^ ^ {Cuu + Cdd + Ccc) - ^Css and the charm 
c = Ccc- The Young diagrams ^ show that the states 
in the [20] multiplet are antisymmetric with respect to 
permutation of two electrons and symmetric in the [36] 
multiplet. Notice that the states of the two [4]-multiplets 
describe the combination of an electron in one of the pos- 
sible d, u, s, c states with an eh symmetric neutral state 
involving the other electron Qo = 5 {dd + uu + ss + ccj ; 
the two electrons are in a mixed state, i.e., do not have 
any speciflc permutation symmetry with respect to one 
another. The two flavor [4]-multiplets are combined with 
their conjugate coordinate parts to form an antisymmetric 
state as dictated by Fermi statistics. The state Qq is the 

flavorless unitary singlet -, which carries no SU4 quan- 
tum numbers (all zeros). This state is also the spin- and 
isospin-singlet with S — and S = 0. 

The structure and the multiplicities of the obtained 
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Fig. 3: (color online) Bound X~ states and their multiplicities 
for (a) u — —2 + e and (b) v = 2 + e. Only the seed k — 
states are shown. Energy zero corresponds to the CR mode at 
huuc- The absorption spectrum of the bright state at Jz — 1/2 
with energy — 0.09i?o is shown in Fig. [T] 

single particle oscillator quantum numbers 7n < 15 are 
used and diagonalize the Hamiltonian matrices numeri- 
cally. The latter typically include about 800 basis states, 
depending on the particular value of Jz. The form of 
the Coulomb matrix elements and other numerical details 
were presented in [T5]. We performed calculations for tri- 
ons with no spin or pseudospin flips for the fllling factors 
i^ — jj, zL e with /i — —2,-1,..., 2. Stable bound X~ 
states were found only for /i = ±2. Our results are valid 
in the presence (as indicated in Fig. [5]) as well as the ab- 
sence of spin Shujs and valley Shuiy LL splittings, provided 
Shujs^v ^ Eq. In very strong magnetic fields B > 20T, the 
n — LL splitting may become of the order of O.lE'o [S]; 
most probably due to the spontaneous symmetry breaking 
[TJ[2]. If this occurs, the SU4 multiplets wfll split, which 
requires additional group theoretical analysis. 

The SU4 symmetry allows one to determine, in a very 
direct and simple manner, the selection rule for trion 
photocreation hio + e'' -^ {eeh). Indeed, the inter- 
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action of electrons with light of frequency uj and cir- 
cular polarizations a is described by the Hamiltonian 
SH± = Al^Q^^j, where t± = t^ ± iTy are the isospin 
Pauli matrices acting in the space of the graphene sublat- 

C —iujt 

tices, A = ^"^rj^ and £, the electric field amplitude 

[TS1I211I25]. When expressed in the SU4 terms, 5H± as- 
sumes the form A (dd + uu + ss + ccj ^^ Qq, i.e., is pro- 
portional to the SU4 singlet Qq. This means that the 
photon is fiavorless in the SU4 sense. From this and 
considerations of the dipole transition matrix elements 
{{eeh)\Qo\e~), it immediately follows that the states in 
the [20] and [36] multiplets are all dark while the states in 
the [4]-multiplet may be bright, provided they have proper 
orbital quantum numbers. Let us discuss the latter. 

The optical selection rules for the orbital quantum num- 
bers are Afc = and AJ^ = ±1 in the a^ polariza- 
tions. The former refiects the negligible linear momen- 
tum of a photon in the dipole approximation and the 
fact that k is associated with translations {i.e., is the 
discrete analog of linear momentum for charged states 
in a magnetic field). Let us consider X~ photocreation 
huj„± + e^jj^ -^ {eeh)k j^-k', here and in what follows Jz is 
that of the seed state. The electron in the initial state in 
the n-th LL has jz = \n\—m—-;^, and Afc = reads in this 
case k = m. Therefore, only the X^ states from the SU4 
[4]-multiplet having Jz — \n\ + ^ for tT+ and Jz — \n\ — ^ 
for a^ can be photocreated. We see that the combination 
of the SU4 and orbital selection rules imposes stringent 
limitations on optical transitions. 

As a result, among all of the bound X~ states found 
near the charge neutrality point, the only bright states are 
at J^ = — 2 + e with the binding energy ^ O.IEq {Jz = |, 
active in (t+, see Fig. [21[a)). Depending on the broad- 
ening, this may give rise to a separate absorption peak 
below the CR or may show up as a low-energy CR shoul- 
der. This is shown in Fig. [1] where the spectra are broad- 
ened ^ 0.05Eo, which amounts to 2meV a.t B — 10 T. 
We conclude that observation of the bright X^ feature is 
almost feasible in graphene even for the current broaden- 
ing/disorder experimental parameters such as mobilities of 
'^ 17000 cm^/Vs [SI ITSII^ . We may compare the binding 
energy to that found for trions in the conventional two- 
dimensional electron gas (2DEG). This is also ~ O.lii^o 
[20] ■ but Eq is smaller than it is in graphene, due to a 
larger value of the dielectric constant e for the 2DEG. 

The oscillator strength of the bright X^ with Jz — ^ 
needs some clarification. Indeed, the hole-like collective 
excitations at 1/ = — 2-|-e originate in the lower cone (inset 
in Fig. [T]). As such, they are optically inactive in the 
electron-like a'^ polarization. This holds, however, only 
in the lowest order in the Coulomb interactions [18]. In 
fact, the resonant ei ^ -'^00-1 processes described by the 
Hamiltonian Hl\ and indicated in Fig. Ulead to "dressing" 
of the optical response. As a result, the oscillator strength 
of the bright X~ becomes comparable to that of the CR 
mode per electron. Thus, the intensity of the X~ peak 




Fig. 4: Photon polarization operator 7r(a;). The box shows the 
resonant excitation of the e-h pair with the deexcitation of the 
electron to n = LL. Because of the multiple Coulomb inter- 
actions (wiggly vertical lines) the intermediate state becomes 
a bound trion Xj^_j, with two electrons in the n = LL and 
a hole in the n — —1 LL, leading to a singularity in n{ijj). 



relative to the main CR is proportional to the number of 
excess electrons e p7H5D] . Notice that because of the e-h 
symmetry, there exists the X~^ counterpart at — i/ = 2 — e, 
which is bright in a~ polarization. 

The dark bound X^ and X^ states may in principle be 
observed by tunneling spectroscopy, which is a sensitive 
tool for probing discrete states in the spectrum. This can 
be pursued, e.g., in tunneling experiments involving gate- 
tunable graphene quantum dots [3T]. Also, the various 
symmetry-breaking terms such as external fields, disorder 
and lattice defects, ripples, and deformations [TJ[2] may 
partially lift the limitations following from the orbital and 
SU4 selection rules. This can make some of the dark states 
"gray" and detectable in the absorption spectra as well as 
possibly in photoluminescence |32) . 

In conclusion, we considered Coulomb correlations in 
graphene with a low density of excess electrons or holes in 
partially filled LL's. We have shown that bound negative 
X~ and positive X+ trions exist in the spectra. With the 
SU4 symmetry of graphene these are analogous to quark- 
quark-antiquark complexes. They are different from their 
counterparts in the conventional 2DEG systems such as 
GaAs [?7H5n] . We expect that the predicted additional 
absorption peak due to the bright X~ below the CR mode 
should undergo an evolution with filling factor e of a par- 
tially filled level. This collective excitation can serve as 
an optical probe as its behavior should reflect correlations 
in the electron system [351|33]- This evolution would be 
especially interesting to follow and study in the vicinity of 
e = i and other FQHE fractions. 

After completion of this work we became aware of pub- 
lication 34 , where the observation of a plasmaron excita- 
tion, a zero field analogue of the X+, has been reported. 
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